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ABSTRACT 


A computer program for simulating the performance of 
mee LOC]B Surface Effect Ship is used to study the longi- 
tudinal motions of the ship under various wave conditions. 
fn investigation into the effect that waves have on bow 
seal, stern seal, and plenum presures is conducted. The 
relationship between the different pressures and their 
associated requirements of input air from the fans is 
Pauoaed., it 25 concluded that the computer program is 
Himited to a certain range of speeds. The concept of the 
wep Capvuring air to reduce drag and increase its speed 
is questionable due to the vapid replenishment Od S266 


mequired to keep the ship riding on its bubble:‘of air. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


aerodynamics 

bow seal 

fourier analysis 

deroedvyaemec force in x direction 
bow seal force in x direction 


BuoUeke—-Gype force #0 x direction due to dracr 
between bubble and wave surface 


Pugider force in x direction 

svern seal force in x eure ion 
Sandee slorce sia, x direction 

mass moment of inertia about x-axis 
mass moment of inertia about xz-axis 
mass moment of inertia about y-axis 
mass moment of inertia about Z-axis 
imieval condit lem 

integral 

Poke moment 

length of craft (72 feet) 

mass 

(ess (16 iGialss fslilie Weibllejemls 

pitch moment 

rotational velocity about x-axis 
pressure in the bubble 

roll angle 


propellar 





PSF 


PSI=wW 


SAM 
SIDWL 
STNSL 
p 


a 
T 


THETA=6 


U 


<< =~ SS < 


Pounds per square Loot 

yaw angle 

meraciomal velocity about y-axis 
air flow rate in 

Ae SielOw abe Ou 

rotational velocity about z-axis 
right hand side 

Shear and moment 

sidewall 

Swieimiel secill 

standard atmospheric reference 
arbitrary steady-state time 

igre Lik feyayeull= 

fonei tudinal velocity 

lateral velocity 

vertical velocity 
foumezonbalwdrescvance in direction of motion 


horizontal distance perpendicular to direction 
of motion 


vertical distance 





IT. @NTRODUCTION 


Conventional monohull displacement ships are limited to 
a maximum speed of approximately 50 knots. This is because 
the drag on conventional displacement ships increase exponen- 
tially with speed, resulting in a limit as to maximum 
meenom ical propulsion plant that can be installed. 

Many new designs of ships have been proposed over the 
past few years. These include the "slender body" ship, 
@ecamaran, hydrofoils, semisubmerged ships, air cushion 
vehicles (ACV), and surface effects ships (SES). Ali of 
these types of ships or vehicles have certain advantages 
and disadvantages. 

One of the most promising designs in overcoming the 50 
knot speed limit is the surface effect ship. The surface 
effect ships are usually classified as to either the Captured 
Air Bubble (CAB) type or the Hovercraft (GEM-Ground Effects 

Machine) type. 

| ~The Hovercraft degen is size limited due to the 
inefficient use of forced air to keep the craft above the 
water surface. No attempt is made to "capture" the air 
Mmaer the craft although the design of flexible skirts 
around the craft have been improved in recent years to 
retard the leakage rate of the forced air. 

The Captured Air Bubble design uses rigid sidewalls and 


flexible seals (bow and stern) to "capture" the forced air 





fmene plenum, thereby reducing the size of the machinery 
memuered TO generate the input air flow. The drag created 
byeuhe CAB craft is reduced considerably when the craft is 
operating in the "on bubble" mode. This in turn allows for 
Mmemeaver speed than conventional monohull ships for the 
~emermamount of thrust. The CAB is not size limited as is 
imees case of the Hovercraft. 

mewdate three CAB type vehicles are operational. 

a. The XR~-3 craft was built by David Taylor Model 
Basin (now designated Naval Ships Research and Development 
Center [NSRDC]). It has a displacement of approximately 
2.5 tons. The XR-3 is presently located at the Naval 
‘Postgraduate School (NPS), Monterey, California. 

b. The SES 100~A craft was built by Aerojet-General 
memeporavion. It has a displacement of approximately 100 
tons and is located at Tacoma, Washington. 

Gee Ene soho. FO0—B erait was built by Bell Aerospace 
Company with headquarters at Michoud Station, Louisiana. 

Ite has a displacement of approximately 100 tons and is 

presently located at Lake Pontchartrain, Louisiana. To 
date this craft has been run at speeds in excess of 70 

Meets (Ref. 1]. 

To better understand the motions and loads of the SES 
meo-B CAB, the Surface Effect Ships Project Office in 
meehington, D.C. issued a contract to Oceanics, Inc. of 
New York to write a computer program and gather certain 


masie Gata on the SES J00-B. The computer program was 
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completed in 19/71 and a copy of the program was installed 
meeene W.R. Church Computer Center at Monterey in October, 
wo 2 . 

mes Invesuipation will concern itself only with pre- 
femeane results obtained in exercising the computer program 
developed for the SES 100~-B CAB. 

imeamestvuer GO Doctor G.J. Thaler, NPS, Monterey, from 
fame. Davis, SESPO, (Ref. 2], a list of ante Studies of 
Mmaemono LOO-B craft, using the computer program developed 
mammpoceanics, was suggested. This author selected to study 
m@emvarliation of pitch, heave, and pressures in regular 
waves (ahead and astern) by varying ship speed, wave length, 


and wave height as suggested in paragraph 2a, [Ref. 2]. 


IESE 





II. NATURE OF THE PROBLEM 


fe LEADING PARTICULARS 
mimes itollowing is a list of the leading particulars of 
Mgemwors LOO-B craft. 


ie, Dimensions 


Length, Overall Or a) tie 
Beam, Overall Son Giles 
Hemp. “Overall 1 | 220. See 
Hull Length {iO Ear 
Cushion Area 1907 ie. 
Bubble Pressure 92.8 PSF 
Distance Centerline to Sidewall aye Sy ee 
Plenum Length at Water Surface OSs ton 
Plenum Width at Water Surface Slee eon ce 
Plenum Average Height Gay, ft 
Froude Number 0.6 

Keel to Deck TRS io, hee 

2. Weights | 

Empty Weight . Hea. O71 lbs 
Gross Weight 209,999 lbs (105 tons) 


beeeeciver Of Gravity 


Longitudinal Feet Forward of Transom So Gas 
Transverse Feet to Starboard Or 0 wel. 
Vertical Feet above Keel 109. Tt: 


eZ 








1 L o 





PEATEEN 
iv | | 


ae 
_ 


EROowLS .On 


Three (3) Pratt and Whitney FT 12A-6 Marine Gas 


Turbines. 4500 HP (each) Rated at 59°F. 


ProOpenmlers: 


Two (2) Supercavitating CRP Marine Screws, 


Diameter 


Feet Forward of Transom 


Feet from Centerline to Propeller Center 


necwenbove Keel 

Fans 

Number of Bow Seal Fans 
Speed of Bow Seal Fans 
Number of Cushion Fans 
Speed of Cushion Fans 
Number of Stern Seal Fans 
Speed: of Stern Seal Fans 
Maximum Total Fan Power 
hucaers 

Number of Rudders 


Feet Forward of Transom 


Feet from Centerline to Rudder Centroid 


Feet Above Keel 
Rudder Area (one) 
Appendages (fins) 
Number of Appendages 
Appendage span 


Appendage Chord 


3 


Ba Ot es 
-1.0 ft. 
1S cal Se ae 

WON) Adee 

i 

L7Oe RPM 
2 

eu RPM 
2 

1870 RPM 

TI Sy nals 

2 

Oe bant. 
1S, 1c eae, 
eS. ft. 

eS 1s 

2 

oe. One 

So) ia 





Distance of Appendage Centroid Forward 

Gio lieeansom te 0 tate. 
Distance from Centerline to Centroid 

of Appendage OS to. 
Distance of Appendage Centroid Above 

nee =a) ft 
Appendage Cant Angle (positive cant 


angle is inboard) | 30.0 deg. 


B. COORDINATE SYSTEM 

Pmeoeooy moving in a fluid can move in all six degrees of 
freedom of motion -- translation along each of three 
orthogonal axes and rotation about each of the three axes. 
The SES 100-B craft has a plane of symmetry (as have most 
ships in existence), which is defined by the ship's center- 
mane and a line perpendicular to the ship's deck. The axis 
system describing the SES 100-B craft is assumed to trans- 
mauve With the craft and only to allow moration of the axes 
about the vertical axis. With the inertial frame axes being 


a right handed cartesian coordinate system, represented by 


X >» ¥ 


E, and Zs with Z POSstuIve down, and an origin selected 


AG 
Bema reference point on the undisturbed free surface, the 
Z-axes are identical. The ship axis system has its origin 

at the center of gravity of the ship with the x-axis positive 
toward the bow, y-axis positive to starboard, and z-axis 


positive downward. The relationship between the inertial 


reference frame and craft axis system is based upon the 


14 





foemies of pitch (8), roll (¢), and yaw (¥) angles. This 


relationship can then be expressed [Ref. 3] 


xX cos ¥Y Sin O xX 


O 
q, = |-sin ¥ cos ¥ 0 Ne (1) 
i 0 O 1 vg 

0 : aul 


C. WAVE AMPLITUDE 
the computer program does not allow for an instantaneous 
build-up of wave amplitude. Instead an exponential rise 


in wave amplitude is used based on the formula 


- time 


Ly AO2% (2) 


wave amplitude = 1.0 -e 
A one foot amplitude will require approximately 4.0 
seconds to build up to 95 percent of maximum amplitude and 
approximately 6.0 seconds to build up to 99 percent of 
maximum amplitude. The reason for this exponential build-up 
or wave amplitude is that the integration routine used in 
Meewcomputer program cannot calculate accurately such an 


instantaneous jump. 


D. COMPUTER SUBPROGRAMS 
ihe computer program developed for the SES 100-B craft 
uses various subprograms for calculating data. These 


Subroutines are [Ref. 4]: 


lbs, 





1. Main Program 


The main program contains the logic for inter- 
connecting the various subroutines. It compares the running 
mome Of time with the finish time, calculates the next 
meme Of time for printing after each print time, and 
Provides some output data. It contains the logic for 
integration of the trajectory (x,y,¥) by the trapezoidal 
male. 

PireUcrwo LIMe, Us Vo W, P,Q, RR; 2, Phi, Theta, 
fer, Si, Sideslip angle, rudder angle, and thrust. 

2. subroutine AEROD 

Subroutine AEROD calculates the aerodynamic 
Sereces and moments on the craft. 

Output is aerodynamic forces and moments. 

Se oubroutine BOWSL 

suoeoutame BOWSL scalculates the forces and moments 
mame cO the bow seal and the leakage flow rate associated 
feet aly Gaps under the various bow seal stations. 

Output is the gap at the various stations, wetted 
dengths at the various stations, and forces and moments due 
mo che bow seal. 

4. Subroutine COLFIL 

Subroutine COLFIL gives a tabular summary of 
selected variables in the vertical and/or lateral plane. 
A Looe of eight variables can be selected for each 


summary. 
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5. oubroutine DMINV 
Subroutine DMINV is a matrix inversion package used 
mer inverting the mass matrix. The method used is the 
meanoard Gauss—Jordan technique. The determinant of the 
miorax 1S also calculated. 
There Is no output from this subroutine. 
6. Subroutine FAN 
subroutine FAN is used to calculate one iano LO 
mes plenum from the fans for a given pressure differential 
@eross the fans. 
Output is flow rates and pressure differentials 
memoas the bow seal, plenum, and stern seals. 
foe Fimvetion FG1 
mime purpose of FGl is to perform the task of 
evaluating a tabular function y = f(x). The procedure 
used is linear interpolation between data points. 
iieromicemo GuLput data from function FGl. 
8. Subroutine FORIT 
subroutine FORIT is used to obtain a Fourier analysis 
Paea periodically tabulated function. 
Memo ie OULPUL Gata from subroutine FORLT. 
9. Subroutine INCON 
Peeroumnes INCON contains the logic for the reading 
Seeall input data, the initialization of variables, and 


mor Che initiation of new cases. 


ie/ 





Stem ronaecompiete listing of all initial condi- 
meens. The initial conditions will always be printed at 
ime Start of each new case. 

A sample data deck for establishing initial 
@eno2tions is contained Pero penaLx 8. 

me. subroutine INTGRL 

puwlbroutcine INIGRL is used to integrate a system of 
first order ordinary differential equations. It uses a 
weerable time step technique based on the Runge-Kutta-—Merson 
meee ichm. The program will stop the calculation if the 
time step becomes smaller than 107° SEC. 

OuGputs will occur when the calculated error is 
meeauer than the tolerance set on that particular integra- 
meee Output format contains time, integrator step size 
(DELT), integrator number (J), error value, and tolerance 
Walue. 

ie DUDRMeUt Ime PROP 

subroutine PROP calculates the forces and moments 
on the fesuuaue, GO tae propulsion system. The port engine 
is always assumed to be the engine which fails if an engine 
Out is designatéd. 

Output is the propulsion forces and moments. 

me. oubroutine RHS 

Subroutine RHS contains the FORTRAN expressions for 

the right hand side of the system of first order differential 


equations. 
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CiepNume@nsicus Of plenum pressure, fan power re- 
G@uered, fan flow rate, leakage flow rates for bow, sidewall, 
and stern seals, plenum area and volume, total forces and 
ements, bow and stern accelerations, and an array of 
Values of the right hand side of the differential equations. 
| mon 6Goubroutine RUDDER 

morourtime HUDDER calculates the forces and moments 
@ememe rudder as well as the rudder motions. 

Output is the forces and moments on the rudder. 

14. Subroutine SAM 

svigenvane SAM performs all the processing and 
calculations involving determination of shears and moment. 

This subroutine was not required to be used in 
@eeainang Che data contained in this report. 

ime oubroutine SIDEWL 

SupZoutsne SIDEWL calculates the forces and moments 
Gue to the sidewalls as well as the leakage flow rates 
Pereectaved with any gaps under the sidewalls. The subroutine 
also contains tice lCewentor the forces and moments due to 
the appendages (fins). 

Output as the gap and immersion depth at designated 
Stations along port and starboard sidewalls and also the 
forces and moments due to the sidewalls. 

mS. Subroutine STNSL 
subroutine STNSL calculates the forces and moments 


acting on the craft due to the stern seal as well as the 


IMS 





leakage flow rates arising from any gaps which open under 
mee seal. 

Cuepub 1S gap and stern seal wetted lengths at 
Mmemerols Stavions across the stern seal in addition to forces 
Pm@aemoments arising from the stern seal. 

meee functions Tl and Te 

Poncumens lL and T2 are two function subroutines 
meread tO calculate various trigonometric relations used 
by subroutine WAVES. 

Mitten omniG OuLpuL data from Functions Tl and Te. 

18. Subroutine WAVES 

Spubroutine WAVES calculates the wave forces and 
Memenus acting on the craft. It also generates the wave 
@mptitudes at various stations around the seals and sidewalls, 
meer l as bubble volume lost due to wave elevation. The 
Meeesovacve for irregular seas is computed by adding together 
a series of regular waves with an appropriate distribution 
Samplitude and frequency. 

Output is the wave elevations relative to calm water 
eeevune Various stations around the bow seal, stern seal, 
Pere and starboard sidewalls. Also the wave elevation at 
the center of gravity, plenum volume lost due to waves, 


and wave forces and moments are listed. 


me) EQUATIONS 
The ten (10) first order ordinary differential equations 


feerved (Ref. 3] for use by a digital computer are: 
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: =1 home nea Glue 
wos mz eS 
p = !o sos ee So 
se, (7) 
ps aoe UN (8) 
wz (9) 
Pe } (10) 
S| = 6 Gris} 
m = p.(Q,, - Q 4) (12) 


In addition three (3) auxiliary equations [Ref. 3] are 


used. These are 


te 
= f u dat (13) 
O 
‘E 
Ve—ieny dt Cy 
O 
t 
Y= fr dt (ILS) 
O 


since the variables x, y, and ¥ defined in Equations 


(13) - (15) are expected to have (in general) slow variation 


ea 





mmmdo Often extend to large values, Simpson's Rule 
Meennique is used for integration in the Main Program 
eeui0n. 

Equations (3) - (12) use the Runge-Kutta-Merson method 
memmi@erration. A tolerance level is selected for each of 
Mm@oe eECGQuations. A truncation error, €, is calculated. If 
Sac tement in € > tolerance level, then the trial time 
meemeoize previously selected is halved and the computations 
are repeated for that time step. If all elements of 
fueely 16 tolerance, the step size is then doubled in the 
memo incvepration step. Otherwise the computations proceed 


meme une Same step size for the next calculations. 


F. NUMERICAL INTEGRATION TECHNIQUE 

mere are many integration schemes to choose from for 
a general system simulation where speed and accuracy is of 
paramount MORE ance . The Runge-Kutta-Merson method was 


selected [Ref. 3]. : 
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ice SCOMEUTER PROGRAM PROCEDURE 


A. PROGRAM FAMILIARITY 

The first priority in using the computer program devel- 
Secanby Oceanics, Inc. was to become familiar with the 
memos suDbrouLiInes and to determine the limitations of 
Miemprogram. Very little data is available at present on 
m@empertormance of the SES 100-B; either through computer 
@umeariOn Or actual craft tests. Since the computer 
meeeram was installed at NPS, Monterey in October, 19/72, 
Mmmeomauchor was one of the first persons to use the program 
installed in the NPS 360/67. As a result, this author was 
fren a growing lad; first learning to crawl, then walk, 
weemranally to run as more and more knowledge was obtained 


meearding the computer program. 


Pee MULTIPLE CASES : 

When making multiple cases (runs) using one data deck, 
all control statements not specified for that ersboe aol lehe 
case will carry the values of the specified data from the 


previous control statement over into the next case. 


fee GIMETATIONS 
The limitations that were determined [Ref. 4] prior to 
use of the computer program were 
1. Wave Components 
The maximum number of wave components that can be 


Heeaq is ten (10). 
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ee Wind 
Neomerovlsvone's made in the computer program for 
iemoanes wind direction and wind speed. 
Em Speed 
Tiewmiianum epeed that can be used is the hump speed 
(16.32 kts) which is where the craft goes into the "on 
bubble" mode. 
4. Water Contact with Top of Plenum 
MemprovasitoOn 1S made in the computer program for 
mememrauving the effect of water striking the top of the 
mem. ne program does print the time, immersion depth, 
and location along sidewall where this water contact does 
eeeur. 
5. Draft 
ii range Of drafts must be selected so that the top 
of the plenum does not sink below the water surface. For 
the SES 100-B, a maximum draft of 6.17 feet is allowed. 
6. Fan Flow Rates and Fan Power 
Fan flow rates and the associated fan power required 
TO produce the fan flow rates are considered to be instan- 
Caneous values. No delay times are used in the computer 
program when fan power and fan flow rates are changed in 


value. 


eee ENETIAL CONDITIONS 
Tg eeerst daca that must sbe obtained is the initial 


conditions for various speeds under calm water conditions. 
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No waves are used so control statement 01102 (see Appendix 
Eemse omitted. In addition the speed is a constant input 
aemecne thrust is allowed to vary by placing a one (1.0) in 
columns 26-35 of control statement 00105 (see Appendix B). 
To insure that calm water steady-state conditions are 
reached, the problem time should run to at least 40 seconds. 
The initial steady-state variables that are required as 
Moe data are draft, pitch angle, plenum pressure, and 
meet. This data can be obtained by placing a one (1) in 
emer, 35, and 55 on the IBM card following control 
statement 00102 (see Appendix B). This data could be made 
mere Ool the variables printed under the vertical and lateral 
plane summaries of control statement 00104 (see Appendix B). 
ieeminitial condition values obtained are shown in Table I. 
dhe ratio of problem time to computer time was found to be 
tees for establishing initial conditions. 

mee check to see if initial conditions are correct, a 
zero (0.0) is placed in blank 26-35, or leave blank, of 
control statement 00105. A blank is Coulee reGgmuanZ ero 1) 
Miemcompurer problem. This establishes the value of thrust 
as initial condition and allows the speed to vary. The 
wemuwes Of the initial condition variables should not change 


with time if the initial conditions are selected correctly. 


me TOLERANCE VALUES OF INTEGRATORS 
The next item that must be established is the tolerance 


value to be assigned to each integrator. The recommended 
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TABLE I 
INITIAL CONDITIONS 


(AS DETERMINED UNDER CALM WATER CONDITION) 


ereeBD U DRArT Je aed e ss Thiol «ONE ENG 
Cais ) CeayASEC) ena} (IN) (DEG ) CL BIS) 
io. 5 PTE 1 AN 17.04 0.09 Sl FAO rs 
20.0 pee78  waho6 16.87 0.12 6490.6 
Zo). 0 ee 22 We eal eh Si ‘Ohya ss. 5 S| ye 
50. 0 50.67 Inn sia lowed Oss 4882.7 
pe .5 54,89 1.334 ioe Ol 0.19 4832.4 
33.0 5 oes 13 0 5. 96 Ore 20 4832.3 
Bo. O sisi lal Wg SEES 1s fae G2 |. 4860.2 
40.0 67.56 lec Wee 15526 Owes 5091.4 
45.0 WG. OL eee 5 Ne Sets: 0.26 5490.2 
50.0 84.45 ilo G 13.99 0.28 S002 <2 
p>. 0 92.90 i103 eeu O32 Li ~6589. 3 
60.0 DOs oe 033 2240 ORS Ss 1 220d 
65.0 109.81 0.956 alee Oe Sys Te 6 
Ao. 0 ILilis (22 Oe ovat 10.49 0.36 eS2in. 
fa. O LAG sarde 0.786 9.432 Oy as Ore Fai 
20.0 Gita 0.694 8.33 0.39 Ceo fF 
85.0 ie 6 1 0.598 7.176 0.39 = 10254.3 
90.0 Heo =O 0.498 5.98 0.40 10697.3 


Note: Plenum pressure of 92.80 PSF is used as initial 


condition for all speeds. 


26 





procedure [Ref. 4] was to set all values to the order of 10 
18. 1072, Miciecati sume Valine of one tolerance, rerun the 
mememem, and compare output data values. If the results do 
Mmoeemcnange drastically, then double the tolerance again. As 
@eamees are noted in the integrator outputs, hold the toler- 
ances for those integrators fixed and vary other tolerances 
Seerreche largest tolerances possible are obtained. 

mes anpucs and output of each integrator are shown in 
eee. o1nce simulation runs were made into or with waves 
mamemesoveady course, it was found that integrators number 
three (3) and ten (10) caused the greatest change in output 
M@eeeeior a factor of ten (10) change in their tolerance. 

Meee selected values of tolerances for each integrator 
mmaeewere used in all computer simulation runs are listed 


iieoable IT. 


Atsueybslel, TAR IE 


TOLERANCE VALUES FOR INTEGRATORS 


INTEGRATOR NUMBER TOLERANCE VALUE 

iL 0001 

2 owe 

3 .000001 

\ .0000001 

5 .0001 

6 .000001 

7 .000000001 
8 .000000001 
G SOOO 
10 .0001 


ay 


» 
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Runge-Kutta-Merson Method 


RKM 
SR 


Simpson's Rule 


Peiuewiamtnipues and Outputs of Integrators 
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mie Values Of the tolerances listed in Table II gave 
output data that had less than one (1) percent difference 
in values when the tolerance values were decreased by a 
factor of ten (10). If integrator number three (3) or ten 
(10) was increased in value by a factor of ten (10), the 


output data would change by as much as six (6) percent. 


Seen itTTAL PARAMETERS SELECTED 
The desired range of parameters initially selected to 
Obtain the data requested by paragraph 2a [Ref. 2], are 
imesved below. 
1. Speed 
pecea ranges Of 30 knots to 90 knots were initially 
selected in 10 knot increments to be investigated. 
2. Wave Direction 
Waves directly ahead and astern were selected to be 
mevestigated. 
3. Wave Length 
Wave lengths from one-half lamda (1/2iA) to five 
lamda (5X) were te ially selected to be investigated. A 
is equal to the length of the SES 100-B which is 72.0 feet. 
4, Wave Amplitude 
Wave amplitudes from 1 foot to 6 feet were initially 


selected to be investigated. 


G. CHANGING WEIGHT 
Symmetry of weight distribution is assumed between port 


and starboard side of the SES 100-B when using the computer 
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program. A nodal method is used to specify the weight and 
Meee Ot gravity of each node. Control statement 00202 is 
Msea with the nodal method. Appendix C lists the various 
mmo Values for the nodal distribution of weight and center 
re ecravity. 

ieeeeca1tional weight is to be added to the craft, the 
femme ot the weight to be added should be halved since 
merry 1S used. The value used as input will only repre- 
memo ial amount added to the starboard side. The center 
mmeravyiry specified will be the same distance to starboard 
evo port. 

Senero! Statement OO020] uses the overall weight, center 
Meee ravicty, and mass movement of inertia about x-axis, y- 
fees, 2-axXis, and xz-axis. 

myen though control statement 00201 meee only one IBM 
Seeen as COmpared to the numerous cards used with control 
statement 00202 (see Appendix C), it was found from numerous 
memmeuver runs that the computer time decreased by approxi- 


mecely 10 percent when control statement 00202 was used. 


H. PROBLEMS ENCOUNTERED 
1. Minimum Speed 
The first problem encountered in using the computer 
program for the SES 100-B was that the longitudinal velocity 
meyedid not reach a steady-state value when waves were used 
at a speed of 30 knots. Varying wave lengths from 36 


feet ( 4X ) to 360 feet (5A) and wave amplitudes from 
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Meeueeieet to 1.0 foot at a speed of 30 Knots caused the 
Memereudinal velocity to decrease exponentially as Tees te 
Maot! hump speed was reached. This result showed that the 
mmmemowdinal acceleration did not go to zero in order for 
the longitudinal velocity (U) to reach a steady-state 
@@ndition. 

fn investigation into the various forces in the 
Meemreccvion was conducted. Table IIIT is a list of the calm 
water force components in the x-direction for various 
Speeds. Force component FXPWAV decays exponentially as 
e156 with an increase in speed. The other force compo- 
nents (FXBS, FXSS, FXSW, FXRUD, FXAED) increase as us with 
Sumeecrease in speed. FXP and Thrust are equal and opposite 
fmrerue and represent the sum of all the forces in the 
Meee ction. A plot of FXPWAV and the other forces in the 
x-direction versus speed (Fig. 2) shows that the cross-over 
occurs at approximately 34.0 knots. 

PMoaelivesuieation into the minimum critical speed 
Oe was sonceene (see Appendix A) with the result that 
U was Eepprowamavely equal to 31.0 knots. This result 


oc 
1.56 and uc 


agrees fairly well with the cross-over of the e 
forces. Based on these results, the minimum speed used to 
obtain the output data was established at 35 knots. 
e. Maximum Wave Amplitude | 
A 2 foot amplitude (height) wave is the same as a 


4 foot high wave measured peak to trough. Wave amplitudes 


of 6 and 8 feet were tried. In this particular case, the 
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Sempucer program would stop in approximately / to ll seconds 
Memmeoolem time Gue to the size of the incremental integrator 
step size becoming smaller than 1 x 107° seconds. Varying 
mmemcolerance values assigned to each integrator would not 
Seeminate this condition. The result is that the Runge- 
ieee erson method of integration is restricted to a certain 
fee, Ol wave lengths and wave amplitudes in order that the 
error generated, when compared with tolerance values, be such 
Zeoeete incremental step size of integration would be halved, 
but not halved until the incremental step size was less than 
1 x 107° seconds. Some limits of speed, wave length, and 
mecnamplitude are given in Table IV, for both ahead and 
Peeett Waves, where the computer program stopped due to the 
@eeeeemetival integrator step size becoming less than 
ip 4 1076 Tcelcacme Uneecvery case, this condition occurred as 
eee emlc Of integrator number 3. 
BemeWater Contact with Pienum 

fiemconpuuer progarm Will print out the time, immer- 
Sion depth, and Teese aioy along the sidewall whenever the 
amplitude of the wave would strike the top of the plenum. 
No provision is made in the computer program to account for 
the additional forces created by the wave striking the top 
of the plenum. It was assumed that immersion depths of only 
meveral inches above the top of the plenum would have very 
little effect on the performance of the craft. Since the 
top of the plenum is 6.17 feet above the keel, immersion 


depths equal to or less than 6.6 feet were considered as 
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mevane negligible effect. Data obtained where the immersion 
depth was greater than 6.6 feet was considered invalid. 
Some limits of speed, wave length, and wave amplitude are 
meen in Table IV for both ahead and astern waves. 
4. Negative Fan Horsepower 

imum Lane morsepower 1S calculated based on the 
emeurm>o OL fan inflow rate into the bow seal, plenum, and 
stern seal. Under certain wave conditions, a negative fan 
mrnow race was indicated for either the bow seal, plenum, 
Siovern seal. Normally the total fan inflow rate would 
meee irom 3,000 cubic feet per second to 10,000 cubic feet 
Memesecond. Several different computer output runs indi- 
Gaved that the fan inflow rates to the bow seal would vary 
Memeween 1,000 cubic feet per second to -100 cubic feet per 
second. ‘This small amount of negative fan inflow rate, when 
moemearead Fo the large net positive value of total fan inflow 
_emereewas considered not sufficient to alter the output data. 
When the net total fan inflow rate was negative, which 
normally petuires a negative fan inflow to all three loca- 
tions, eee tavive cToOval fan horsepower resulted. Since the 
fans can not be put in reverse (pumping air out), this 
condition was considered severe enough to cause the output 
@eea tO be invalid. Some limits of speed, wave length, and 
Wave amplitude, for both ahead and astern waves, are listed 


in Table IV where negative fan horsepower resulted. 





Je04 9°9 UU 

JasieyT yAdeq UoOTSA9WUT 
adamodesdoy uey esaztqeson 
190g 9°9 UeYY 

JlasizeTyT uadagq uoTtsiawwy 
IamModasJoYH uey sATAeDSON 
eal! UB 

Ssoy7 ozts deiac s.oeigequy 
04 enp paddoaras wedisorg 
geet UeBYY4 

ssoeq eztc deis Jo4eida4jul 
04 enp peddois weidgodg 
_OT4T UbUY 

SsSoy azts Perec Saaseas ae 
O04 enp peddoris wedsodg 
TaMOdaesiJoy uez sATAeBsan 
JeamModesIoy uet saAT yegon 


TJahnodssioy ueg_ sat Aegon 


(reg ret 


(idagdeT 


(IssieT 


(AagIeT 


(TIgIeT 


pug) 


pue) 


pue) 


pug) 


pue) 


O° 
Cie 


O° Tf 


mes 
O°SL 


0° Or 
0°09 
O°SE 


O°SE 


O0°SE 
0° Or 
0°09 
0°SL 


en tn pti i a i x = hr A 


uoSsPedoYy 


SHALANVEVd PIAVSN-NON 
aTaVia 


~- 


(bel 
epn4yttduy 
JACM 


0°9e UL94Sy (SSeT pue) 
eo) ere ULBISY (IodIe@T pue) 
O°80t UI94SY (SSOT pue) 
0°80T pesouy (IagieT pue) 
OF 9e. peeuy (IeZIeT pue) 
Oe peouy (laddie—T pue) 
O° O09 pesuy (losieT pue) 
O°9E peouy (IadIeT pue) 
O°9e peouy (Ies1eT pue) 
O°8Ol peouy (TadJe—T pue) 
(14) (SLY) 
U4sUIsT HO bald o.rG poads 
OACIM oACM TeRIFUL 


36 





5. Problem Time Versus Computer Time 

The ratio of problem time to computer time for a 
1 foot amplitude wave was found to be approximately 1:40. 
[iiereesing wave length and speed increased this ratio very 
Meeenely (1:42). This ratio holds good for both ahead and 
Peeemiewaves, With just a slight decrease in the ratio for 
astern waves (1:38). Increasing the wave amplitude from 
ieome feet would change the ratio from 1:40 to 1:60. 
eereasing the wave amplitude from 1 to 3 feet would change 
Geemratvio from 1:40 to 1:70. 

6. Reducing Computer Time 

HG was found that using the calm water initial 
conditions would require approximately 140 to 160 seconds 
memune value of U to reach a steady-state condition. With 
famaooitional 10 to 20 seconds added for printing out steady- 
Meme values, a total of 2.90 to 3.90 hours of computer time 
Wee required for each combination of speed, wave direction, 
Weave length, and wave amplitude. 

To reduce the aes ,- time required, a run of 40 
Peconds problem time was used with the initial calm water 
Pemaioions aS input values. The value of U obtained for 
each run was plotted between 0 and 40 seconds and an esti- 
mate of steady-state speed (U) was made. This steady-state 
speed (U) could be estimated because the speed (U) would 
decrease exponentially with time from its calm water value 
as waves were encountered. The new estimated steady-state 


speed (U) was then used in control statement 01201 as a new 
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meee ral Speed, and the program was run for an additional 
emeerereonas. This resulted in a reduction of 1.0 to 1.5 
Melee Ol Computer time for each parameter of speed, wave 
@irec>ion, wave length, and wave amplitude used. 
fee Cucput Data Limitation 

The IBM-360/67 computer at the W.R. Church Computer 
Pemcer has a built-in limitation as to the maximum number 
Saeeees Of Output data that can be obtained.’ A rule of 
mmmb used was that a maximum of 100 time intervals of 
femeecOUulLG be printed out if all the print switches listed 


under control statement 00102 were ON (indicated by a 1.0). 
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IV. PRESENTATION OF DATA 


Pen AG VERSUS SPEED 

wetee ore ctrar versus Speed under calm water conditions 
memeiown in Fig. 3. The results show that a minimum value 
(saddle point) in the curve occurs at a speed between 32.5 
Seems.) KNOLS. Comparing Fig. 3 with the results reported 
Mamenei. 5, close agreement does exist between the value of 
hump speed and the value of drag occurring at hump speed. 
Peeare 5 Shows that the minumum value of drag occurs ata 
[emer speed and that this particular value of drag is less 
mare conat reported in Ref. 5. 

mpeeure 3 indicates that once hump speed is reached and 
mmiemcrart goes "on the bubble", it would accelerate from 
16.5 knots to approximately 68 knots under calm water 
@emaivions if the value of thrust was held constant. 
Reference 5 Imei oaeed Dia oeenererentewould accelerate trom 
hump speed to a speed of approximately 40 knots if the 


[ue OL thrust was held constant. 


B. STEADY-STATE SPEED 

Bach particular wave condition resulted in a different 
steady-state speed being obtained. Table V lists the steady- 
state speeds that were obtained from various wave directions, 
wave lengths, and wave amplitudes. Although the steady-state 
Value of speed is only valid for the particular wave condi- 


tion listed, the values are useful in estimating other 
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Meemay-svacve speeds under wave conditions that are not 
fsved. This will reduce computer usage time if the proce- 


dure listed in Paragraph III H-6 is followed. 


C. BOW AND SIDEWALL LEAKAGE 

Mmcer certain wave conditions, air leakage would occur 
from the bow and/or sidewalls. Of course this is in addi- 
MmuemeeOo Tne stern leakage that occurs Mader fe sied 
Memeeoi1ons. To better understand when leakage from the bow 
and/or sidewall did occur, Table VI oe Che steady-state 
Seem speed, wave direction, wave length, and wave amplitude 


when this leakage occurred. 


Peer nASE RELATIONS OF VARIABLES 

To help understand how the simulated motions of the 
craft would respond to varying wave conditions, a phase 
comparison of selected variables with respect to wave motion 
Meseconducted. It was found that changing the .speed of the 
meee Qld not alter the phase relation between the selected 
variables and penile wave used under both ahead and 
astern Wavemorreceions. im addition, changing the wave 
amplitude did not alter the phase relation between the 
selected variables and the particular wave used under both 
ahead and astern wave directions. <A representative initial 
condition speed of 60 knots was selected to show the phase 
relation between the different variables with respect to 


various wave lengths. 
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MmeenvVariables that have been selected for study, based 
Memeeet. 2c are: Draft, Plenum Pressure, Bow Acceleration, 
Memeer Of Gravity Acceleration, Bow Seal and Sidewall 
Leakage Rates (where applicable), Stern Seal Leakage Rate, 
ips Fan Flow Rates (bow, plenum, stern), and Total Fan 
Horsepower. 

ij Ahead Waves 

Wewemilenetis of 36 feet (1/2 A), 72 feet (A), ana 
ieemereet (3/2 X) are shown in Figs. 4 through 24. 

Figures 4, 5, and 6 show that the response of the 
Memver Of gravity acceleration leads the wave for a wave- 
length of 36 feet and lags the wave for wavelengths of 72 
and 108 feet. The smallest magnitude of center of gravity 
@eeeweracvion occurs for the {2 foot wavelength. Figures 
7, 8, and 9 show that the bow acceleration is in phase with 
Mmemecencer of gravity acceleration and its magnitude is 
memeer. his is to be expected for a rigid type ship. 

iiewaverate  Gramt response is shown in Figs. 10, ll, 
ogemle. No distinct phase relationship exists since the 
atic s of draft are averaged along the length of the side- 
wall. Plenum pressure is the primary determinant in estab- 
jishing average draft values. The higher the plenum pressure, 
the smaller the average draft value that will occur. The 
T2 foot wave length has the smallest draft with the corre- 
sponding largest plenum pressure (see Fig. 14). 

Figures 13, 14, and 15 show the response of plenum 


Pressure to different wave lengths. Pressure in the plenum 


a 





mamomiinarily based on the plenum volume and the difference 
Meeween the air flow in and out of the plenum. As the 
volume is decreased (increase in draft), the plenum pressure 
mmeeeimerease. Comparing plenum pressure with average 
feo. you can see that the two are in phase. In addition, 
Meemeoressures into the bow and stern seals are in phase 

mei vie plenum pressure. 

Figures 16, 17, and 18 show the stern seal leakage 
Mmeeeroccurs under dadifferent wave lengths. The value of 
stern seal leakage is based upon a calm water baseline 
[eaeece rate of 6417 cubic feet per second at a plenum 
Mressure of 92.8 pounds per square foot. Since it is not 
known precisely where all the leakage does occur in the 
Crait, it is assumed that all of the baseline leakage does 
mummerOouc the stern seal. The values of stern seal leakage, 
Under wave conditions, are based upon the one iatelalyo) Iie enh 
Venturi flow. For the same area of stern seal gap, the 
amount of stern seal leakage will vary proportionally with 
che plenum Puesoure, SAS autvesult, the stern seal leakage 
Maeve 1S in phase with the plenum pressure. No bow or 
sidewall leakage occurred at the wave lengths shown. 

ome sinpub fan flow rates to the bow seal, stern 
wemeeanGa plenum are shown in Figs. 19, 20, and 21. The 
fan flow rates are based on a function of pressure. Since 
Mgemitan Characteristics are nonlinear, a series of linear 
curve-fits for the pressure versus input fan flow rate is 


used. The slope of the pressure versus input fan flow rate 
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Mmeemeees at difierent locations on the plot, with the result 
fee vne fan flow rate causes a more rapid increase in fan 
ieew rates at lower pressures. Basically the fan flow inputs 
meemmeto OL phase with the pressures. In Fig. 19 and el, 
Meeecsve tian flow inputs to the bow seal are not shown. 
meenmmaximum negative value that did occur was approximately 
meme pounads per square foot. 

Figures 22, 23, and 24 show the total fan power 
required to produce the necessary input of air flow. Of 
interest here is to note the double frequency of oscillation 
That occurs. The total fan power is based upon a summation 
Sarat) POwer required for the bow seal, plenum, and stern 
Peete in turn each individual fan power is based upon the 
Mieessure times the fan flow input. Since pressures and fan 
Meret nputs are out of phase, the maximum power requirements 
Wall occur when the pressures are the largest and also when 
meee) ft 10w inputs are the largest. This results in the 
eemmere frequency of oscillation of the total fan power 
Requirement. 

2. Astern Waves 

Wave lengths of 72 feet (A) and 108 feet (3/2 d) 
are shown in Figs. 25 through 40. Wave lengths of 36 feet 
are not shown for astern waves since they follow approximately 
the same sinusoidal phase relation established for the 36 
foot wavelength of ahead waves. 

Figures 25 and 26 show the relationship of the 


center of gravity acceleration to the wave. The maximum 
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acceleration lags the wave by a small amount. Figures e7 
and 28 show that the bow acceleration is in phase with the 
center of gravity acceleration, as is to be expected for 
meet d type ship. 

iemaveratenmcaral, response 1S Shown in Fig. 29 and 
30. No phase comparison between draft and wave can be made 
Since this is an average value. The value of draft is 
Memiarily dependent on the value of plenum pressure. 

Figures 31 and 32 show the response of plenum 
Pressure to different wave eee ns Basm@ally tic plenum 
pressure follows the average draft response, but Fig. 32 
emews Now the plenum pressure can deviate from a sinusoidal 
mereror. Both bow seal and stern seal pressures are in 
meee With the plenum pressure. 

Figures 33 and 34 show the stern seal leakage for 
G@@iaerent wave lengths. The stern seal leakage rates are 
in phase with the plenum pressure response as a result of 
the same reasoning used for ahead waves. 

Figures 35 and, 36 show the bow and sidewall leakage 
Petes. No sidewall leakage occurs at the 108 foot wave length 
(Fig. 36). Leakage will only occur when a gap is created 
between the bow seal and water surface and/or the sidewall 
and water surface. The amount of leakage rate is dependent 
Upon the area of the gap created and the pressure of the 
eenum. The principle of Venturi flow applies to the leak- 
age rate. When sidewall leakage did occur (Fig. 35), this 


moeaee Occurred at the forward part of the sidewall. 
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The input fan flow rates are shown in Figs. 37 and 
38. As with ahead waves, the input fan flow rates are out 
mame se With the pressures. 

Figures 39 and 40 show the total fan power required. 
@_meemi@orice the double frequency of oscillation that occurs. 
This occurs because of the same reasoning used for explaining 


the ahead wave case. 


E. AVERAGE AND FLUCTUATION VALUES OF VARIABLES 

mpyeragse and fluctuation values of selected variables 
Memmevar lous speeds, wave directions, wave lengths, and wave 
amplitudes are shown in Figs. 41 through 80. The variables 
meemeare Shown in the different figures are: Pitch, Center 
eames ty Acceleration, Bow Acceleration, Draft, Bow and 
mee ocal Pressures, Plenum Pressures, Stern Seal Leakage 
Rate, Bow Seal Input Fan Flow Rate, Plenum Input Fan Flow 
faee, obern Seal Input Fan Flow Rate, and Total Fan Power. 
meemeructuations, and not the average values, of Bow Acceler- 
ation and Center of Gravaty Acceleravion are shown Since 
their average values are approximately zero. 

The average values were calculated by adding the values 
Stee number of data points (20 to 50) along a full wave 
Emeounter cycle and dividing by the number of data points. 
Either one or more full wave cycles were used depending on 
the number of data points available in each wave cycle. 
This summation of values was necessary because of the 


non-sinusoidal nature of certain variables. 
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The values of fluctuations used were based on taking 
the difference between the maximum and minimum values of 
mee Variables. 

The dots shown in Figs. 41 through 80 represent the 
Gata points. No interpolation was used between the data 
points. 

The values of wave length and wave amplitude are denoted 
Meee Ligures at the right hand side of the last data point 
Meeq. A 36 foot wave length and 1 foot wave amplitude is 
denoted as 36 xl. 

Figures 41 through 60 are for ahead waves (directly on 
the bow) and Figs. 61 through 80 are for astern waves 
(directly astern). 

A Berscmssion of the results obtained in Figs. 41 


meoueh 80 is contained in later sections. 
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Vee DISCUSSION 


PemeechrITURE AIR BUBBLE CONCEPT 

The Capture Air Bubble (CAB) ship is designed to run 
along the water surface on a cushion (bubble) of air. The 
Semeepe 1S to "capture" the air supplied by fans in the 
plenum by using rigid sidewalls and flexible bow and stern 
meee Lhe air that is captured in the plenum area will 
@eemree the draft of the craft and in turn will reduce the 
meee Ol drag. This allows a higher speed at a given 
Smee OL thrust for the CAB ship than a conventional 
Memon! craft of the same gross weight. 

imeone computer Simulation of the 1OO-B SES, an initial 
base leakage of 6417 cubic feet per second was assumed to 
Peer OuL the stern seal under calm water conditions at a 
plenum pressure of 92.8 PSF. The plenum volume, under calm 
feieeeeeOnditlons, varied between 9,200 and 10,200 cubic feet. 
The variation in plenum volume is due to the speed Oi the 
Meet. At slower speeds Piewdiad aids eeveaver Causing 42 
decrease in plenum volume. 

Nommeacctcer What the speed of the craft is, the calm water 
results indicated that the volume of the plenum was replenished 
Bpproximately every 1.5 to 2.0 seconds. A leakage of 6417 
Spore feet per second may not be considered as extremely 
excessive, but the fact that the plenum is replenished 
fer. 5 to 2.0 seconds appears to be in contradiction to 


the theory of the "captured" air bubble. 
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Peeures 51 and /1 show that increasing speed and increa- 
sing the wave amplitude will reduce the average rate of 
stern seal leakage. The average plenum pressures and the 
mead Stern seal pressures also decrease with an increase 
in speed and wave amplitude as shown in Figs. 47, 49, 67, 
and 69. This agrees with the computer program since leakage 
Mmeeemas dependent upon the area of the gap created and the 
meeeoclaved pressures. However, the average input of air 
Mmeomune DOW seal, stern seal, and plenum also decreases with 
an increase in speed and wave amplitude as shown in Figs. 
eee i, (3, (5, and 7/7. This is contrary to the computer 
program where a decrease in pressure causes an increase 
iieiiputc air by the fans and vice versa. The average total 
memeepower required also decreases with an increase in speed 
and wave amplitude as shown in Figs. 59 and 79. 

The computer program assumes an initial gap at the stern 
me of 25.52 square feet. An additional gap is calculated 
due to wave conditions. This additional gap will occur only 
when the wave surface is approximately below the level of 
the Peel. Hine total gap at the stern is then based on the 
Miitial gap and any gap due to the wave. It was found that 
mie area of the gap due to the wave had little or no effect 
Meche total area of the gap at the stern. This was pri- 
marily due to the small amplitude waves used in the study. 

f@omethe results obtained of the input of air by the 
means and the total fan power required, it shows that the 


more severe wave conditions require less input of air and 


des 





MeescpOower. The volume of the plenum was still being 


replenished by air about every 2.0 seconds. 


B. WAVE CONDITIONS 

Wave lengths of 72 feet (A), both for ahead and astern 
Memeo, Nad the smallest effect on the performance of the 
100-B SES. The amount of decrease in speed from the initial 
speed input to the steady-state speed was the:smallest for 
ment OOL Wave length. In addition, the center of gravity 
Seeeweravion and bow acceleration were less for the 72 
foot wave length than either the 36 foot or 108 foot wave 
Memethn. Generally the average stern seal leakage rate, 
meme tal power, bow seal input of air, plenum input of 
_meeocern seal input of air, and bow seai, stern seal, 
maemo lLenum pressures were higher. 

meerourph it is not known why this result did occur, it 
appears that a wave length approximately equal to 1G Oe 
memeod Of the craft gives the best performance of the 


@eait in regard to minimum motion. 


ee 





VI. CONCLUSIONS 


A. PERFORMANCE 

Mme periormance of the SES 100-B, as determined by 
computer simulation, did not vary appreciably between ahead 
(directly on the bow) waves and astern (directly astern) 
waves under the same wave conditions. Steady-state speeds 
Maeeemeareer, average plenum pressures were slightly higher, 
the fluctuation of pitch angle was larger, and the average 
draft was slightly less for astern waves as compared to 
ma@eaa waves. The other performance variables shown in 
Paes. 41 through 80 were of similar values. 

Meer similar wave conditions, the computer program 
@eewea that the craft could proceed through larger amplitude 
Meeern waves than ahead waves. Table IV showed, for 
example, that a wave length of 36 feet and wave amplitude 
of 2 feet from astern could be’ run at a speed of 60 knots 
(vice only 50 knots if the wave was from ahead) before a 
Merpacive fan power developed. The crairt Commas run ava 
speed of 40 knots in an astern wave of wave length 36 feet 
and wave amplitude of 4 feet before an immersion depth 
greater than 6.6 feet occurred. At the same speed, a 
negative fan power condition developed for an ahead wave 
of- wave length 36 feet and wave amplitude of 3 feet. Wave 
lengths and wave amplitudes in this range are associated 


Merely a state 2 to state 3 sea condition. Therefore, it 
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fears that the CAB concept, and particularly the SES 100-B, 
femeeoe Linited to calm or moderate sea conditions. 

When bow and sidewall leakage did occur, generally the 
meme Ol leakage was larger at the slower speeds for a 
Biven wave condition. Table VI shows when bow and sidewall 
[eeeeeree Cid occur. Increasing the wave length resulted in 
larger leakage rates out the bow and Sidewall. Under a 
given speed and wave condition, the magnitude of bow and/or 
Sidewall leakage was larger for an astern wave than an 
eifeaa wave. 

As speed was increased, in general the average values 
mmeeeoech and draft increased with a corresponding decrease 
in their fluctuation values. The fluctuations of the center 
Sueeravivy acceleration and bow acceleraticn also increased. 

As speed was increased, in general the average values of 
Meetea!, Stern seal, and plenum pressures decreased with 
Smameicrease in their fluctuation values. The average 
fees Of Stern seal leakage also decreased, in general, 
feed COrresponding increase in fluctuation as speed was 
Miereased. The average inputs of air by the fans to the 
mevesecal, plenum, and stern seal decreased with a corres- 
Memories increase in fluctuation as speed was increased. 

The same trend also held for the average and fluctuation 


values of total fan power. 


B. RECOMMENDATIONS 
The following recommendations of additional areas of 


Study are made on the basis of the results obtained during 
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irs Study. 
im Leakage Rates 
The resuits of the computer program showed that the 
fees Deing replenished in the plenum every 1.5 to 2.0 
feeemds. Hrrom the CAB concept, this appears to be an 
excessive amount of air input required. It is recommended 
that the results of air leakage be compared against those 
Seeome OFS 1O00-B to see if this is actually ae case. 
pee Air Inputs by Fans 
The results of the computer program showed that an 
average decrease in air input by the fans and an average 
G@ecrease in the amount of fan power occurred when a more 
Meeere Wave condition was encountered. This does not 
appear feasible from the CAB concept. It is recommended 
maeeea Study of the modeling of the air inputs by the fans 
and leakage rates be conducted to resolve this apparent 
discrepancy. : 
3. Stern Seal Leakage 
The amount of stern seal leakage is based upon the 
Mees) pap (both initial gap and gap due to waves) that is 
ereated at the stern seal. The amount of stern seal leakage 
memonen dependent on the total gap and the plenum pressure 
mimeo exists. It is recommended that the computer program 
be revised so that the amount of stern seal leakage is 
listed in two parts: leakage due to the initial gap and 
any additional leakage due to a gap being created by the 


waves. 
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4, Steady-State Speed 


A minimum initial speed of 35 knots was used based 
Samene results obtained in Fig. 2 and Appendix A. No 
Bempucer runs were made to obtain the minimum critical 
speed between 30 and 35 knots under different wave conditions, 
Meemeanead and astern. It is recommended that an investiga- 
Pomeincvo the minimum critical speed that can be used in 
meee computer program be nr to determine what effect 
Miererent wave conditions will have on the values of this 
minimum critical speed. 

meee validity of Computer Data 

beetle or no performance data on the SES 100-B is 
[atlable here at the Naval Postgraduate School. To check 
meemvalii0ity of the computer data obtained in this study, 
mecomparison of the performance data of the SES 100-B 
Soomro be made under similar sea conditions as used in 
Maes scudy. A complete evaluation of the computer program 
Meeeeie selected variables that were studied could then be 
made. 

5. Additional Data 

ine range of wave conditions was very limited in 
Memeo seudy. Although basic trends could be developed, it 
ieerecommended that additional computer runs be made using 
Simple waves from abeam and then using complex waves composed 
mame Or more different waves. Realistic sea conditions 
mrould be modeled by the use of multiple simple waves as 


mopuUt data. 
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APPENDIX A 


CALCULATION OF CRITICAL SPEED 


From Newton's Second Law it is known that 


f= ma Cad) 


meee longitudinal force, re is equal to 


Pee hruse — drag (252) 


Ponsicering thrust as being constant, then we need only 
memooncerned with those terms involving u to find the 


miemptcal value of speed. Since u is composed of those 


2 -1.56 


terms that are u”~ and u dependent, Eq. (3.1) and 


Pome .3-c) can be written as 


ae = =i ue + Kk, yrt+ 56) + Thrust (3.3) 


Linearizing Eq. (3.3) by letting u = uy + 6u gives 


‘ 2 2 
~nm §6 U = k, (uy + 2udu + 6u) 


i eae 6u,-1.56 
+ ky Uy OR oe ) 


O 


= Thrust (3.4) 


Peolying linear approximation, since we are interested 


mea steady-state critical speed for Uys and by speciltying 


PS 





that °° = — Eq. (3.4) can be written as 


6u 


many = 3/2 
(1 + us) 


~3/2 : 
- ecese 


-méu = (2k, u,) é6u + kou 


Letting x = ou » vhen 
O 


ae ~ Were 
; a | 


2 (1 + x) 


hoping that 


I. 1 
= oF ee ee > ae (3.7) 
VE Oe RN 2 See 


(1 + x 


i 
(l + x)(l + x)72 


and that the value of x, which is equal to 2%, will be 


O 
much less than 1, Eq. (3.7) can then be written as 


a ee. (3.8) 
(1 + x)(1 + x) Cecbee etal) 


Thus Eq. (3.5) can be written as 


~5/2 


-md6u = (2k, us? su - kU, éu 


[(2k,u,) - ead Su (3.9) 
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For 6u = 0, the critical speed, Use? iomicom fo. (5.9 ) 


2k.u = kyu “5 fe (3.10) 


Poeevine for eae gives 
ee 
areal = em Coa) 


The values of k. and kj~ are determined from the calm 


oD 2 
feeeerei1orces listed in Table III for 30 knots. 


FXPWAV = 5.898 x 102 


my 
li 


Woraaleror cose meee 667 x 10° 


te 
an 
iH 


eee ae ek and Kj into feet per second gives 


| C 


Be O98 xX 10° 
a GonGn).-/* 


_ 3.867 x 10° 
1. (50.67)° 


mere 50.607 feet/second = 30.0 knots. 


Pouvane Ba. (3.11) for Lee gives 


Uno = 51.5 feet per second = 31.0 knots. 


ILS, 
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